Late-time entropy production due to the decay of domain walls 
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It is shown that late-time decay of domain walls can dilute unwanted relics such as moduli, if 
the universe was dominated by frustrated domain walls with tension a = (1 ~ 100 TeV) 3 . Since 
energy density of the frustrated domain walls decreases as slow as the inverse of the scale factor, 
an overclosure limit on the axion decay constant f a is also considerably relaxed. In fact f a can be 
as large as the Planck scale, which may enable us to naturally implement the QCD axion in the 
string scheme. Furthermore, in contrast to thermal inflation models, the Affleck-Dine baryogenesis 
can generate enough asymmetry to explain the present baryon abundance, even in the presence of 
late-time entropy production. 

PACS numbers: 98.80.Cq, 11.27.+d 



Particle physics beyond the standard model predicts 
a number of new particles, some of which have long life- 
times and decay at cosmological time scales. If such long- 
lived particles are substantially produced in the early 
universe, they may give crucial effects on cosmology and 
spoil success of the standard big-bang model. 

One well-known example of such dangerous relic par- 
ticles is gravitino in supergravity theories. Gravitinos 
are produced during reheating after inflation and the 
abundance of the produced gravitino is proportional to 
the reheating temperature. The gravitino with mass 
~ 0.1 — 10 TeV decays during or after the big-bang nu- 
cleosynthesis (BBN) and may destroy the light elements 
synthesized in BBN. Thus, in order to keep success of 
BBN the reheating temperature should be sufficiently low 
(for recent works, see pi y> L» H ) ■ 

Another class of dangerous relics are light scalar 
fields called moduli which appear in superstring theo- 
ries 0, 0, ■ The moduli fields are expected to acquire 
masses of the order of the gravitino mass from nonper- 
turbative effects of the SUSY breaking and have long life- 
times since they have only the gravitationally suppressed 
interactions. During inflation the field value of moduli 
generally deviates from the vacuum value due to extra 
SUSY breaking by the cosmic density. The deviation is 
of the order of the Planck scale M G (= 2.4 x 10 18 GeV) 
since no other mass scale exists. Then, when the Hub- 
ble becomes less than the modulus mass, the modulus 
field starts to osciilate with amplitude ~ Mg. Because 
the modulus denity is comparable to the cosmic density 
at onset of the oscillation, it soon dominates the den- 
sity of the universe and causes the serious cosmological 
problem (moduli problem). Since the abundance of the 
moduli is independent of the reheating temperature, the 
moduli problem is much more serious than the gravitino 
problem. To solve the moduli problem (as well as the 
gravitino problem) it is necessary to dilute the oscillat- 
ing moduli by huge entropy production. So far the most 
successful model to produce such large entropy is thermal 



inflation proposed by Lyth and Stewart |a li| ■ The ther- 
mal inflation is mini-inflation with the number of e-folds 
~ 10 which takes place at the weak scale. A comprehen- 
sive study of the thermal inflation was done in Ref. 01 
and it was shown that for modulus mass between 10 eV 
and 10 TeV the thermal inflation solves the cosmological 
moduli problem. However, the thermal inflation also di- 
lute the baryon number of the universe and hence we need 
an efficient baryogenesis mechanism. In Refs. 0, it 
was pointed out that the Affleck-Dine mechanism does 
that job. But later it was shown that Q-ball formation 
obstructs the baryogenesis [l^ . 

In this letter we propose an alternative model to di- 
lute the dangerous cosmological relics. In our model do- 
main walls, which are usually considered as cosmological 
disaster, play an important role. The domain walls are 
produced when some discrete symmetry is spontaneously 
broken. The domain wall network in the universe be- 
comes very complicated if there are many discrete vacua. 
The density of such domain wall network decreases very 
slowly and quickly dominates the universe. In extreme 
case, the domain walls become frustrated and their den- 
sity pow decreases as pdw ^ o, ( a: the scale factor) 
which is slower than the scaling evolution pow °c a~ 3 / 2 
in matter-dominated universe. Furthermore, if the scalar 
potential has a tiny term which explicitly breaks the dis- 
crete symmetry, the vacua separated by domain walls 
have slightly different energies and hence the domain wall 
network is no longer stable and decays producing large 
entropy. It will be shown that the entropy production by 
the decay of domain wall is enough to dilute the moduli 
density. Moreover, because the density of domain wall 
decreases much more slowly than the matter density, the 
cosmic axion density is also diluted. As a result, the ax- 
ion decay constant f a as large as Mq is allowed, which 
might enable us to implement the axion into the string 
framework. In addition, the Affleck-Dine baryogenesis is 
compatible with our model as we will see below. There- 
fore, the domain walls provide more efficient and consis- 
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tent dilution mechanism than the thermal inflation. 

First we present our scenario that the decay of domain 
walls generates large entropy to dilute unwanted relics 
such as moduli and gravitinos, and derive constraints on 
tension of the domain walls. In contrast to the other can- 
didates for late-time entropy production, domain walls 
have an advantage that the energy density decreases rel- 
atively slowly, which ensures huge entropy production. 
To be specific, an equation-of-state w equals to —2/3 if 
the domain walls are completely frustrated, that is, the 
structure of the domain wall network remains unchanged 
as the universe expands. Throughout this letter, we as- 
sume this is the case. Later we give a model which would 
lead to such frustrated domain walls. 

Let us consider the evolution of the energy density of 
domain walls, which eventually dominate the universe. 
The energy density of the domain walls at the formation 
is estimated as 



PDW,i 



a Hi, 



(1) 



where a is the tension of the wall, Hi the Hubble param- 
eter at the phase transition. We have assumed that there 
are only a few domain walls per one horizon when they 
are formed. Here and in what follows we neglect O(l) nu- 
merical factors. If the domain walls are fully frustrated, 
the energy density falls off as the inverse of the scale fac- 
tor: pow a -1 . For definite discussion, we assume that 
the universe is dominated by either the inflaton or mod- 
ulus fields when the domain walls are formed. Noting 
that the energy density of an oscillating massive scalar 
field decreases as oc a -3 , the domain walls dominate the 
universe when the Hubble parameter becomes equal to 



H eq ~a-HtM G 



(2) 



If the reheating occurs earlier than H = H eq , the 
radiation-dominated epoch might last for a while, which 
makes it easy for the domain walls to dominate the un- 
verse. Thus, for conservative discussion, we assume that 
the universe is dominated by nonrelativistic particles un- 
til the domain walls dominate the universe. 

The decay of domain walls can proceed if there is a 
tiny bias between vacua, Sp. When this energy difference 
becomes comparable to the energy of the domain walls, 
the decay occurs and the universe is reheated. Since the 
universe was dominated by the domain walls in our sce- 
nario, 5p is simply related to the decay temperature as 
Sp ~ T%. In order to have both large enough entropy 
production and the successful BBN, we take the decay 
temperature as low as lOMeV. Note that this value is 
just an exemplified value and that the successful dilution 
is actually realized for a wide range of the decay temper- 
ature, say, T d = lOMeV - lOGeV (see Eqs. ©, Q and 
{7J)). It is the smallness of Sp that enables domain walls 
to be long-lived. Let us note that, when the domain walls 



decay, at least several domain walls must be present in- 
side one horizon, otherwise the old inflation occurs some- 
where leading to unacceptably inhomogeneous universe. 
This requires that the decay temperature Td satisfy the 
following inequality: > aHd, where Hd ~ T^/Mq 
is the Hubble parameter when the domain walls decay. 
Thus the tension is bounded above: 



a < TjM G ~ (100 TeV)- 



lOMeV 



(3) 



On the other hand, the tension should be large enough 
to dominate the universe before the decay: H eq > Hd. 
That is, 



a > 



T%Ml\ 3 



Hi 



(100 GeV)- 



Td 



lOMeV 



Hi 



ITeV 



(4) 

Now let us estimate the the abundance of the mod- 
ulus field after decay of the domain walls. For domain 
walls with the tension satisfying JjJJ and (@J , H eq is much 
smaller than the modulus mass. Since the modulus field 
dominates the universe as soon as it starts oscillating, 
it is reasonable to assume that the universe was domi- 
nated by the modulus field when H = H eq . Then the 
modulus-to-entropy ratio is given by 



Td 



Hd 

Hp, 



o*H ' 



(5) 



where we have used the fact that the Hubble parameter 
falls off as oc a -1 / 2 while the domain walls dominate the 
universe, and Eq. J2J is substituted in the last equation. 
This must satisfy the observational bound: 



< k — = 3.6 x 10 k h GeV, 

s s 



(6) 



where p c is the critical density, K varies from 10 10 to 0.2 
depending on the modulus mass , and h is the Hubble 
constant in units of 100 km/sec/Mpc. Thus the tension 
is further constrained as 



er > K 



■ (500 GeV) 3 



T d V 
10 MeV 



H, 



ITeV 



(7) 



To sum up, the domain walls with the tension a = (1 ~ 
100 TeV) 3 , which decay just before the BBN starts, can 
dilute the moduli to the observationally allowed level. 

Next we present a model that would lead to forma- 
tion of the frustrated domain walls. Let us consider the 
following superpotential: 



N N 



w = ^E E z u*&i + E ( 8 ) 

i j v i 

where Zij and $i are gauge-singlet superfields, and real 
coupling constants A > and e are assumed to satisfy the 
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inequality, A ^> |e| 17]. These interactions respect Sn 
symmetry if $i and Z%j are taken to be the fundamental 
and bi-fundamental representations of Sn permutation 
group. The i?-charges are assigned as i?$ = and Rz = 
2. We assume that always sits at the origin: (Z^) = 
0, which can be realized if Z^ acquires a positive Hubblc- 
induced mass term during relevant epoch. 

Then we obtain the following effective potential for <I>i, 

- V -mlJ2\*i\ 2 - m t/2z2(*i+*?) 

% i 

+A(El^| 2 ) 2 + 4e El^! 4 ' ( 9 ) 

i i 

where m 3 / 2 is the gravitino mass, Vo — 0{m^/X) is cho- 
sen in such a way that the cosmological constant vanishes 
in the true minimum, and we have used A 3> |e|. Here we 
have assumed the negative mass of the order of the weak 
scale at the origin, too ~ h-EW ~ 1 TeV, which is induced 
by the SUSY breaking effects. The third term comes from 
the gravity-mediated SUSY breaking effects. Due to this 
term, the global minima of the potential lie along the real 
axes of $j . In order to study the vacuum structure of this 
potential, let us concentrate on the the real components 
of <£>i by setting Im<I>i = 0. Rewriting the potential in 
terms of the real components 4>% = V^Re^i, we obtain 

V(<f>) * Vo - m 2 J2 + X(£ ^) 2 /4 + e ]T 4, (10) 

i i i 

where the third term in Eq. @ is neglected since we are 
interested in the case of to 3 / 2 ^5 too. The scalar potential 
given by Eq. (|10|) actually agrees with the O(N) model 
studied in Refs. [H 1 1 ll | , The approximate O(N) sym- 
metry originates from the hierarchical couplings, A and 
e. Note that this potential is obtained by disregarding 
Im$i, so it does not necessarily give a right answer, for 
instance, when one estimates the energy density inside 
domain walls, although it is still useful to study the vac- 
uum structure. 

The positions of the global minima of V((f>) depend on 
the sign of e. For e < 0, there are 2N minima: 

^ = (0,...,0,±»i,0,... ) 0), for i = 1 ~ JV (11) 

with vi = mo/yX. On the other hand, if e > 0, the 
potential minima are given by 

(±«2,... > ±t*) (12) 

with V2 = mo / V XN, where arbitrary combination of the 
signs are allowed. In the following, we consider the case 
of e > 0. Then all Z^ fields acquire masses of the or- 
der of 1)2 ~ h-EW m these minima. Also, the number of 
the minima is 2 . Thus, for large N, the vacuum struc- 
ture becomes more complicated, compared to the former 
case. Since there are many types of vacua, pair annihi- 
lation processes of walls are expected to be highly sup- 
pressed, which enables walls to deviate from the scaling 



law. The minima are separated by the potential barrier 
represented by the third term in Eq. ©, if e is larger 
than TOgy 2 /t;| ~ XNm^, 2 /rn^. Then the tension of the 
walls is given by 

(13) 

On the other hand, if e is smaller than to^ 2 /w|, the po- 
tential barrier is the last term in Eq. 10. The tension 
is (72 ~ |e| 5 «2 — where the last equality holds when 
e ~ m 3/2/ w 2- Thus the tension of the domain walls in 
this model can satisfy the bound derived in the previous 
section, although XN J$ 1 might be necessary in the latter 
case. 

We assume that all $j get positive Hubble-induced 
masses therefore sit at the origin until the Hubble param- 
eter becomes comparable to Too- Then those scalar fields 
start rolling down and get settled somwhere on de- 
fined by J2i \®i\ 2 — 2iVuf. When the Hubble parameter 
becomes equal to the gravitino mass, all $^ move to the 
real axes due to the third term in Eq. and fall in 
one of the vacua (|12fl . Then the domain walls with the 
tension u\ or 02 are formed. 

Domain walls must decay before the relevant BBN 
epoch. To this end, we introduce a R-symmetry vio- 
lating interaction that lifts the degeneracy of the vacua: 
W = c' (W) J2i ®1/ m g with (W) = m 3/2M^, leading to 
Va = c 7713/2 Si ^l/Mc+h-.c., where c and c' are coupling 
constants. For simplicity we take c both real and negative 
so that true minimum is given by </>truemin = (^2, ■ • • > ^2)- 
The decay temperature is then expressed by 

lOOMeV v ; V0.1/ VlTeV/ VlTcV/ v ' 

When the bias becomes comparable to the energy of the 
domain walls, the topological defects are no longer stable 
and decay into $-particlcs in the true vacuum. Not to 
spoil the success of the BBN, however, ^-particles should 
decay into standard model degrees of freedom. We as- 
sume that $i weakly couples to some heavy vector-like 
quarks, which enable <i>-particles to radiatively decay into 
the standard model gauge bosons as soon as the domain 
walls decay. Since the i?-parity of <& is even, the de- 
cay processes into the lightest supersymmetric particles 
(LSPs), which may overdose the universe, can be avoided 
if the mass of <E ) is smaller than two times the LSP mass. 
Thus our model of domain walls can naturally satisfy the 
constraints necessary to induce successful dilution. 

As we saw in the above discussions, the domain wall 
can be a viable candidate for late-time entropy produc- 
tion. What differs from the other candidates is that the 
energy density falls off very slowly: pdw oc a -1 . Such a 
distinctive feature can lead to another important cosmo- 
logical consequence: the overclosure limit on the axion 
decay constant can be considerably relaxed. If we do not 
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assume entropy production after axion begins the coher- 
ent oscillation, the axion decay constant is constrained 
as F a < 10 12 GeV not to overdose the universe. This up- 
perlimit is relaxed to 10 15 GeV, if the late-time entropy 
production due to the decay of nonrelativistic particles 
occurs 15]. In the following, we show that the axion de- 
cay constant is further relaxed and can be as large the 
Planck scale if the domain- wall induced entropy produc- 
tion occurs well below the QCD scale. 

The axion starts to oscillate when the Hubble parame- 
ter becomes comparable to the mass: 3H osc ~ m a (T s £ ). 
The axion mass m a depends on the temperature T as [la ] 



rn a (T) 



0.1m a (A QCD /T) 3 - 7 for T>A QCD /^, 
m a for T<Aq CD 7t, 



(15) 

where Aqcd — 0.2GeV, and m a ~ Aq CD /F a is the axion 
mass at the zero temperature. Since the number density 
of the axion decreases as oc a -3 , the energy density of 
the axion when the domain walls decay is 

1 H e 

Pa \ T =T d = ^m a (m a (T osc )F*e 2 )-f 7 (16) 

where we have used H oc aT 1 / 2 when the domain walls 
dominate the universe. 9 ~ 0(1) denotes the initial am- 
plitude of the axion field in the unit of F a . The axion- 
to-entropy ratio is then given by 



Pa 
S 



1.3 x 10^ 



(17) 



e(T osc )5A« QCD M6' 

where £(T) = m a (T) / m a < 1. This must be smaller than 
the present value of the ratio of the critical density to the 
entropy. That is, the axion decay constant should satisfy 

F a < 3.6 x 10 18 GeV£(T osc )i^H 

'qjSy ( Aqcd V ( T d y* 
0.14/ \0.2GeV ) v 10 MeV,/ ' v ' 

Therefore the domain-wall induced entropy production 
opens up the axion window to the Planck scale, if £(T osc ) 
is close to 1. In other words, the axion is a good candidate 
for dark matter if F a ~ Mq. We need to check that T osc 
is smaller than 0.1 GeV so that £(T osc ) is ~ 1. To this 
end, the evolution of the cosmic temperature before the 
decay of domain walls must be specified. If the decay is 
approximated to be an exponential decay with a constant 
decay rate, we obtain 



0.05GeV 



_5lA 

10 MeV J 



0.26 



F a 
M G 



-0.13 



< 0.1 GeV. 
(19) 

10 MeV and 



Therefore, £(T osc ) is close to 1 when Td 
F a ~ M G . 

In this letter we have investigated the dilution of un- 
wanted cosmological relics such as moduli by entropy pro- 
duction of the domain wall network. It has been shown 



that the frustrated domain walls quickly dominate the 
density of the universe and their decay produces entropy 
large enough to dilute the dangerous moduli. We have 
also given a concrete model which leads to frustrated 
network of domain walls. Moreover, we have found that 
the late-time decay of the domain walls greatly relaxes 
the constraint on the axion model and the axion decay 
constant f a as large as the Planck scale is allowed. 

Up to here, we have assumed fully frustrated domain 
walls, however, this assumption can be weaken to some 
extent. Similar arguments show that domain walls whose 
energy evolves as pdw ck a with 7 < 0.2(0.1) for 
k = 0.2(10~ 10 ) works as well. Finally we make a com- 
ment on baryon number generation in the present model. 
The domain-wall decay also dilute pre-existing baryon 
number. Therefore, there should be large baryon asym- 
metry before the entropy production or baryon num- 
ber should be generated after the decay of the domain 
wall. Since the decay temperature is expected to be low 
(~ 10 MeV), it is unlikely to produce the baryon number 
after the decay. Then the most promising candidate for 
baryogenesis is the Affleck-Dine mechanism. In order to 
produce sufficient baryon asymmetry the gravitino mass 
should be relatively small, m 3 / 2 ~ 10 MeV lid , as in the 



gauge-mediated SUSY breaking scenarios [18j. In the 
case of the thermal inflation, the Affleck-Dine mechanism 
does not work due to Q-ball formation [J^. The crucial 
point is that the thermal inflation requires large messen- 
ger scale, which leads to large Q-balls and small baryon 
number in the background universe. However, our model 
allows relatively small messenger scale, since V2 is much 
smaller than the vev of the flaton. Thus the Affleck-Dine 
baryogenesis does work in the present scenario. 
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